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Abstract. We set up the Maxwell's equations and subsequently the classical wave 
equations for the electromagnetic waves which together with the generating source, 
a traveling oscillatory charge of zero rest mass, comprise a particle traveling in the 
force field of an usual conservative potential and an additional frictional force /; these 
further lead to a classical wave equation for the total wave of the particle. At the 
de Broglie wavelength scale and in the classic- velocity limit, the equation decomposes 
into a component equation describing the particle kinetic motion, which for / = 
identifies with the usual linear Schrodinger equation as we showed previously. The 
/-dependent probability density presents generally an observable diffusion current of 
a real diffusion constant; this and the particle's usual quantum diffusion current as a 
whole are under adiabatic condition conserved and obey the Fokker-Planck equation. 
The corresponding extra, /-dependent term in the Hamiltonian operator identifies 
with that obtained by H.-D. Doebner and G.A. Goldin. The friction produces to the 
particle's wave amplitude a damping that can describe well the effect due to a radiation 
(de)polarization field, which is always by-produced by the particle's oscillatory charge 
in a (nonpolar) dielectric medium; such a friction and the resulting observable diffusion 
as intrinsically accompanying the particle motion was strikingly conjectured in the 
Doebner and Goldin original discussion. The radiation depolarization field in a 
dielectric vacuum has two separate significances: it participates to exert on another 
particle an attractive, depolarization radiation force which resembles in overall respects 
Newton's universal gravity as we showed earlier, and it exerts on the particle itself 
an attractive, self depolarization radiation force whose time rate gives directly the 
frictional force /. 
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1. Introduction 

While the usual linear Schrodinger equation has demonstrated to be adequate for 
the common nonrelativistic quantum systems, L. de Broglie suggested [1] in the 50s- 
60s that the quantum mechanical wave equation may be more generally nonlinear. 
Various forms of nonlinear equations have been proposed and investigated subsequently, 
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for a similar concern of internal states of particle as L. de Broglie's or from rather 
different fundamental considerations. Of these, the Doebner-Goldin form of nonlinear 
Schrodinger equation. Doebner-Goldin equation, represents a unique family which H.- 
D. Doebner and G.A. Goldin obtained in [2] by admitting observable diffusion current 
to the probability density of a quantum particle, and subjecting this to the continuity 
equation of a Fokker-Planck type on the basis of a unitary representation of an infinite- 
dimensional Lie algebra of vector fields and group of diffeomorphisms. 

In view that it admits observable diffusion which commonly occurs to a greater 
or lesser degree in all macroscopic processes that are at the microscopic scale in 
majority cases executed by quantum particles, and in view of the physical significance 
associated with the unitary representation of group theory based on which it derives, 
the Doebner-Goldin equation can be anticipated to represent an important prediction 
of certain possible intrinsic processes accompanying quantum systems. What these 
possible processes may be has on the other hand remained as an open question prior 
to the present study. The main question may be formulated as that, what can be such 
an (intrinsic) dissipative process which disrupts not the stationary state of a quantum 
particle as an ordinary heat process would, and which in the meantime manifests itself 
an observable diffusion? Prom a measurement point of view at least such a process is 
viable even to a first degree, since such a process would not cause any detectable effect 
if a measurement is made over the damped probability density current and the damping 
in amplitude does not change with time. A theoretical recognition of such processes 
however would seem unrealistic until recently, in view that the mechanical nature of 
the quantum processes described by the usual Schrodinger equation had remained up 
to interpretation. 

With overall experimental observations as input information the author recently 
proposed [3a-e] an internally electrodynamic (lED) model for simple particles such as 
electron, termed also basic particle formation (BPF) scheme in earlier reports [3a-h, j] 
(with coauthor P.-I. Johansson). The lED model, briefly, states that a simple, single- 
charged particle is constituted of an oscillatory point charge q of a zero rest mass and 
the resulting electromagnetic waves propagating at the speed of light c. In so far as the 
way the mechanism of the model operates, q can be of arbitrary quantity; q is to be 
given as an input (out of two sole input data, the other is the total energy of the charge) 
to yield the actual material particles. [Por examples, of the elementary particles, for the 
isolatable charged one where the multiple- or " neutral-" charged particles are viewed as 
achieved by integration processes n <— p + e + i/g and N <— p + n, clearly then \q\ — e; 
for the, as of today, nonisolatable charged ones, quarks, \q\ = 1/3,2/3.] What form the 
basis of operating mechanism of the model merely are three elementary, experimentally 
firmly corroborated laws regarding electromagnetic waves: the Maxwell's equations in 
respect of wave propagation, the Doppler principle in respect of source motion effect, 
and the Planck energy equation in respect of energy discretization. 

The obvious immediate motivation for proposing the lED particle model was to 
reconcile our understandings of particles with the various puzzling phenomena involving 
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particles to date. For example it was not understood that what are the internal, 
mechanical processes which causes a particle to manifest both as an extensive wave 
and a point object depending on methods of detection? And what are such internal 
mechanical processes which appear to simultaneously also command a particle to emit 
or absorb electromagnetic waves through the exchange of a portion or the whole of its 
own internal energy or inertial mass? What are the origin and nature of mass? In view 
of the uniform presence of charges in all of material particles, of the universality of the 
vacuum as a medium to as far as we know, and of the fact that in ordinary connection 
to charge the only pure waves propagating in this vacuum are electromagnetic waves, it 
is natural to expect that the internal processes of material particles are electromagnetic 
and that the operation of these hold the key to the answers to the various relevant 
puzzles. 

In part as a broad test of the lED model, and in part as an endeavor of the 
understanding of a range of diverse phenomena from a common ground facilitated by it, 
a range of predictions of the fundamental properties of particles and relations have been 
made in terms "first-principles" solutions for the lED partcile internal processes. Here, 
the "first- principles" refer to a minimal set of firmly established physics laws consisting 
mainly the few aforementioned; and the charge and the total (mechanical) energy, — 
corresponding to a characteristic oscillation frequency for a universal vacuum medium, 
are as two sole input data. The achieved predictions [3a-j] include: A particle has a spin 
and relativistic mass[3d-e] apart from the input characteristic charge and total energy; it 
is extensive as the result of its electromagnetic waves being extensive, and when traveling 
freely its waves evolve into a traveling, and in turn a standing, de Broglie phase wave 
between boundaries [3c, e]. As such, the particle obeys the de Broglie relations [3c, e]; its 
traveling de Broglie phase wave will produce constructive interference at integer times 
the de Broglie wavelength upon superposition of its different parts, and it in turn behaves 
like a point object owing to its charge, say, when scattering elastically with another point 
particle. More generally, in arbitrary potential fields under corresponding conditions the 
particle obeys the Schrodinger equation[3d,e] and the Dirac equation[3i]; the particle 
obeys the Einstein mass-energy relation[3c,e], the Galilean-Lorentz transformation[3f], 
and Newton's law of gravitation in attracting another particle[3j], among others. 

It is natural that we extend in this paper the studies to aim to derive the Doebner- 
Goldin nonlinear Schrodinger equation, which we show will result when additionally 
subjecting the Schrodinger particle to a frictional force of the medium with the total 
system subjecting to an adiabatic condition. In part, this derivation provides an 
additional test of the lED basic-particle model. And in part, as the lED model itself 
suggests a natural origin of frictional force to be the radiation depolarization field 
always produced in a dielectric medium along with a particle's internal electromagnetic 
processes, the derivation and the solution will provide a formal elucidation for the 
connection of this with the Doebner-Goldin observable diffusion. 
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2. Particle model with electrodynamical internal processes 

We consider an lED model particle is traveling at a velocity v as its oscillatory charge q 
does, for simplicity in a one-dimensional box along X-axis. The oscillation of the charge 
is associated with a total energy Sg, Sg being smallest at v — 0, Sg\i,=o = tg. 8.g or 
Eg may be endowed e.g. in a pair production in the vacuum. Eg describes the ground 
state and therefore cannot be dissipated or detached from the charge except in a pair 
annihilation. 

The charge will owing to its oscillation generate electromagnetic waves, with the 
radiation electric field E-' and magnetic field B-' (of the j th component) governed by the 
Maxwell's equations given in a medium of dielectric constant k in zero external fields as 

V.E^ = ^, V-B^ = 0, VxB^ = /.f^ + -^, VxE^ = -^. (41) 

Where is the density and j;j the current of the charge q of the particle, assuming no 
other charges and currents present; e = kcq and ~ ^uq) with cq and /iq the permittivity 
and permeability of the vacuum, and c the velocity of light in the medium. (Until 
the specific application oriented denotations in Sec. 6, unless explicitly specified the 
respective media for conveying the particle and for the reference of measurement will 
not be explicitly specified in writing the electromagnetic variables.) Considering only 
regions sufficiently away from the source charge so that — jq — 0) making some 
otherwise standard algebra of the equations (41), and replacing the field variables by a 
more general dimensionless displacement (p^ , with — Aip^ , — E^ /c — A(fi^ /c and 
A a conversion constant, we obtain the corresponding classical wave equation for each 
component electromagnetic wave (f^: 

% = ^^VV^- (42) 

Here, in view of a Doppler effect to result from the source motion (to express explicitly 
later) , we distinguish by the superscript j the component wave generated in the direction 
parallel with the source velocity v, denoted by j — f, and the one in the direction 
antiparallel with -u, denoted by j = |; within walls in stationary state there must 
also simultaneously prevail their refiected components and we may regard these as if 
being generated by a virtual charge (compared to the presently actual charge) which is 
reflected and traveling in — X-direction, denoting by j = virf and j = virj. Apparently, 
the total wave given by the sum of all of the component waves, J2j 9^"' = describes 
the particle. 

Based on the general results of electrodynamics applied here to the particle's 
internal processes as governed by the basic equations (41)-(42), basic properties of 
a given particle can be predicted; in the remainder of this section we outline two 
directly relevant ones of these. The first is the total energy of the wave and accordingly 
the particle. As a general result of classical electrodynamics based on solution to the 
Maxwell's equations, the (41) here, combined with Lorentz force law, an electromagnetic 
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wave (j) transmits at the speed of light c a wave energy and a hnear momentum 
p' = e^/c. In virtue of the stochastic nature of the electromagnetic waves, its 
total dynamical quantities, the total wave energy and linear momentum here, are 
appropriately the geometric means. 



From the underlining laws afore-used, mathematically the amplitudes of E\ B\ (p\ etc, 
and accordingly , p' , E and p are permitted to take on continuous values. 

Following M. Planck's discovery of quantum theory in 1901, it has been additionally 
understood that the amplitudes of these quantities are in nature quantized; the total 
wave energy of an electromagnetic wave of frequency a;/27r is £ = n?iu!, that is, e 
consists in general of n momentum-space quanta, or photons, each of an energy fuj. 
The electromagnetic wave comprising our basic particle, like an electron, positron, etc., 
has, based on experimental indications especially the pair processes, a "single energy 
quantum", n — 1; the Planck energy equation for the total wave of the particle therefore 
is 



This total wave of a single energy quantum here has in a one-dimensional box two 
components, 9?^ and 99*. a situation no different from discussed after (42). Their wave 
frequencies arc Dopplcr displaced to cj^ and c^;* as a result of the source motion (to express 
explicitly below); and similarly as (4), uj'' = e''/h, u;* = e^/h. Further from (49) we have 
uj — y/cJkJ^. For the total wave comprising the particle, s represents therefore the total 
energy of the particle. It has been proven and formally expressed especially through 
quantum electrodynamics that, the Maxwell's equations and naturally the derivative 
classical wave equation continue to hold; and the quantization of the fields and the wave 
energy etc. formally is the result of subjecting the corresponding canonical displacement 
and momentum, corresponding to the u{= aip) and ii here, to the quantum commutation 
relation [-u, u\ = ih. In this generalized framework, clearly the classical solution of a 
continuous amplitude for say e merely is an approximation when n is large. 

The second property is the incrtial mass of the wave and thus the particle. The 
two components of the electromagnetic wave, [E\B^) or rapidly oscillating at a 
geometric mean frequency uj /2t: and wavelength \ — c/ (a;/27r), viewed at some distance 
and ignoring the detail of the oscillations will appear as if being two rigid objects, 
wavetrains, traveling at the speed of light c; the two trains of the component waves 
together make a total wavetrain. In view that its speed of travel, c, is finite as contrasted 
to infinite, the total wavetrain has inevitably a finite inertia mass, denoting this by m. 
This mechanical representation of the total wave, as a rigid "wavetrain", permits us at 
once to express according to Newtonian mechanics the linear momentum of the wavetrain 
to be p = mc. Combining this with the classical electrodynamics result £-|-0 = pc of (49) 
gives the kinetic energy of the wavetrain e — mc^] this is just the Einstein's mass-energy 



£ = P — s/p^'i — y/p^C or £ — pc. 



(49) 



(4) 
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relation (see e.g. [3c,g] for a detailed treatment). This energy and the Planck energy of 
(4) ought to equal of course, thus 

m = huj/c^. {A.l) 

The mass m of the total wavetrain comprising the particle naturally represents the mass 
of the particle, here acquired dynamically through the total motion of the waves of a 
geometric mean frequency ujji-n. m is dependent on the particle velocity and is thus 
relativistic, see further after equation (14) below. 

3. Wave equation of total motion of particle in external fields 

To the particle we now apply a Coulomb force F — — W owing to a conservative scalar 
potential V, and in addition a viscous force /; these give a total applied force F' = F+f. 
We express the / as follows. Suppose out of the total oscillation of the charge, a 
fractional displacement Uq only produces radiation and is defined here to be equal to 
the wave displacement -u = a^, ^ being the dimensionless total wave displacement 
in the field of applied potentials; this in zero applied potential field is the ^0 earlier. 
Making direct analogy to the viscous force of ordinary mechanics, we can write down 
the frictional force opposing the total motion of the particle as / = (^^]^(^^^)"^ ™ 
units of N. That is, in general / is a function of the time rate of the total displacement 
of the charge or alternatively of the resulting wave displacement in the medium; section 
7 will give a concrete representation of such a force. Assuming d{a}^)/dT is small, so to 
a good approximation / = -^-^ = "^^(^ + VV^^), where hi is a constant in units 
of Nms and is real; / is in units of and is generally imaginary (pointed out by D. 
Schuch) for ^ being generally complex. This may rewrite as 

■' ^L^ ^dT dX ' 
where 

iPihvoy,. i/3i;i(vy)V 

2mD 2m|y|2 ' ^ > 

W{= P = duj/dk) is the wave speed of and W* (= -dX/dT = du'/dk') of 
that of the imaginary y*. The expressions in (7) follow firstly from the requirement 
that W is in direct proportion with the velocity Vohs of the current jobs (= '^^obsP) 
of the probability density p {— j'^p) in order to ensure the continuity of current in 
a non-absorbing medium. That is, W = Ivobs — W*, where the imaginary W* is 
subtracted from the generally complex Ivohs- The current jobs = "^obsP of P with a 
uniform translation at velocity v alternatively is according to Fick's law the diffusion 
of a varying p in a viscous medium, jobs = ~-DVp, with D the diffusion constant. So, 
^^obs = ^ = = -fKVy*)^ + y*Vy]; here we have taken D to be as defined in 
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(6) and this will receive a justification later through the role of jobs in (29). The above 
leads explicitly to the first and second expressions in (7) once we put the proportionality 
constant as J = — where (3i is a parameter yet to be determined [by the equation 
(29) below], and the other constants are inserted so that /3i will have the simple solution 
value 1. 

In virtue of the electrodynamic origin of F and inevitably also / which are 
empirically estabhshed for point particles, extending to the extensive lED particle here 
the two applied force and thus their total F' act apparently directly on the point charge. 
We now want to map this F' into a force directly interacting with the E^, or (p^. We 
notice that by its mathematical form equation (42) represents just a classical wave 
equation for the electromagnetic wave ip' , a(p^ therefore a mechanical wave propagated 
in an elastic medium, and F' interacts with aip^ by an effective force -F„jcd acting on 
directly on this apparent medium. On the basis of this direct correspondence, but 
taken as a heuristic means only in this paper (so that we need not to firstly introduce 
with sufficient justifications at any detail the structure of this elastic medium), we 
shall below map the force F' into -Fmed- Now, while F' drives the charge, of a mass 
m of the particle, into an acceleration 9^^/97^, F^^ drives the medium of mass 
(effectively) into acceleration ^^'^med/^T'^ in ^-direction. Supposing the charge and the 
medium oscillate at a fixed phase difference if not in phase, the two accelerations must 
be equal, we thus have 



F' 

med 



m 



{F + f) = -pi 



V 2D fd^ \ 



(8) 



Where OJl<^ = L^pi, with pi the linear mass density of the medium along the wave path 
of a total effective length — JL (^ winds in J loops about the box side L) . 

We below further implement the force F^ed in wave equation (42) similarly using 
the heuristic approach by applying directly Newton's laws to the apparent elastic 
medium. In the apparent medium ^ corresponds to a physical, transverse displacement 
u — — aC J2j V'^ 3.S produced by the disturbance of the charge oscillation, with a a 
conversion constant of length dimension. The deformed elastic medium is consequently 
subject to a tensile force = p^c^, with c the velocity of light at which (/?•' propagates. 
This force and the applied force -F^ed together give the total force acting on the 
particle through directly acting on the apparent elastic medium 



K- Fr- ^Led - Pi 



V 2Dd^ 



(9) 



where the minus sign in front of F^ed represents that this force tends to contract the 
chain. Consider on the linear chain of the medium a segment AL at (X, X + AX) is 
upon deformation tilted from its equilibrium AX an angle 'd{X) and d + A'(9(X + AX); 
assuming ^ is small, will be uniform across the entire wave path L^. The 
transverse {Z-) component force acting on it is AF^^ = F^[sin('i? + A-i?) — sin-i?], with 
[sin(i? + M) - sin^] = [1 + 0{^)\M ~ A^? = V2(ay)Ai?; 0{^) collects the higher 
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order terms and is dropped (this term leads to anharmonicity and not damping, and as 
can be shown this in general leads not to a Doebner-Goldin form of nonlinear term). 
Substituting in the above with (9) for we have 



m y oT mp 



(10) 



Applying Newton's second law to the segment of a mass A9Jt^ = piAL, ~ Pi^X, we 
have piAX^-^^ — AF^^. Placing (10) in it, dividing apiAX, we get the equation of 
motion for per unit length per unit mass density of the elastic chain of medium at X, 
or equivalently the classical wave equation for the total (electromagnetic) wave of the 
particle: 



m )3 oT mp 



In summary, (11) has a basic part = which one will get from summing over 

all j values the wave equations (42) given earlier directly from the Maxwell's equations in 
zero applied potential, and it has an additional part describing the effect of the applied 
force F^ed' derived with the help of the "heuristic elastic medium" approach. 

Concerning the solution of (11), for the present we only consider explicitly the case 
of D = 0. So (11) reduces to ^ = [c^ + ^] V^^g; this being hnear, thus = J2j 
and 



(11) 



9V 



m 



vV' 



(12) 



Supposing also V is a constant, Vc, equation (12) can be immediately solved to consist 
of plane waves, = C cxp[i{k''X — uj^T + 00)], 99* = — C exp[i(/i;*X + cl;*T — ao)]- Where 
y = 'f^K are the Doppler-displaced wavevectors for the wave generated parallel with 
the source velocity v {j = '\) and antiparallel with v {j — J), and u;^ — Q are the 
corresponding angular frequencies, with 7+ = 1/(1 — t'/c), 7* = l/(l-|-f/c). K and 
Q — Kc are the values of and lo^ a,t v — 0, c being the velocity of light as before. 
The explicit superposition of the incident waves ip^, cp^ and their reflected ones 

..J V 



.0. .4- 

,viri ,„vir+ g-^g ^ standing wave (for a systematic representation see [3a-c,e]): 



where 



kd^^{B-K){K-k^)^^Ka, Ka^{^^K- u^V^t^^Q; 



(13) 



(14) 



7 = = 1/ V 1 ~ "^^Z c^- Canceling u between (14) and (A.l) gives further m = 7M, 

with M = hQ/(? the classic- velocity limit {v'^/(? 0) of m, i.e. the rest mass of the 
particle. An explicit inspection of (13) will readily show that k^ is the de Broglie 
wavevector (for an existing elucidation see e.g. in [3c,i]), being its value at the limit 
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v^/c"^ 0. We shall later (see after equation 23) generalize the representation to the 
case where V may be arbitrarily varying in L; until then we shall proceed the following 
discussion for the constant V, Vc- Substitution of ^0 in the total wave equation given 
from the linear sum of the wave equations (12) over all j, i.e. with D = in (11), 
directly gives the expected relativistic energy-momentum relation for the particle[3b], 
which gives an additional check that ^0 of (13) is the correct solution to the total wave 
equation. 

For the solution of wave equation (11) with D finite we shall use the trial function: 

^ = 2^0, where Z = e"^, Q = Qi+iQ2. (15) 

Z represents a damping factor; Qi and Q2 are real variables and are in general functions 
of X, T. We shall restrict ourselves to the case where D is small and accordingly 
\iQ\ << \i{K + kd)X — iu}T\. Under such a condition, for the derivation of a nonlinear 
Schrodinger equation in question below, until the context of equation (33), an explicit 
solution form of Q needs not be known. 

4. Transformation to wave equation for kinetic motion of particle 

At the classic- velocity hmit v'^/c^ — > 0, the total wave function "^^ reduces to (see [3a-c]) 
hm„2/,2^o'y0 = c'e^(^dX-^d^) with f?d(= l^i^T) = ^K^v + which is equivalent 
to the solution for Schrodinger equation for an identical system as described by the 
wave equation (11) in the case oi D = and V = Vc- Therefore, as we noted in [3a-c], 
equation (11) must inevitably have a direct correspondence to the Schrodinger equation, 
and the remaining question mainly then was to identify a physically justifiable procedure 
to transform (11) to a form of the Schrodinger equation. Such a formal procedure was 
elaborated in detail in [3a-c] by a back-substitution of the explicit function "^0 in wave 
equation (11) in the case of D = 0; by use of the Fourier theorem the procedure further 
led to a Schrodinger equation for V arbitrarily varying and also, by a straightforward 
extension, for three dimensions. For the present case of D being in general finite, we 
below similarly first reduce and simplify wave equation (11) at the classic-velocity limit 
v'^/c'^ — s> to a form such that the K- and K^- processes are separable, by means of 
back-substitution of the formal function ^ (15) where the function is explicitly known 
and Q assumed small. 

We first prepare for the separation of the K- and Xd-pro cesses in three aspects, the 
first two being similar as for the case D — [3a,b]: (i) We observe that (11) contains 
the derivative ^ which relates to the acceleration of the particle and, as such, the K- 
and Kd- processes are not separable; but the two processes are separable for the first 
derivative ^ which relates to the total energy (for a detailed analysis see [3b]). This 
suggests us to lower the time derivative one order as |^ ~ ^[(^)^*'^ + ^06**^^1^] — 
^[-iu"^ + 0] = -iiu§. (ii) In the two terms ^V^y and ^^iV^pV^^ in (11), the 
coefficients in front of V^^, being approximately the scale of quadratic thermal velocity 
v"^ or lesser, are relatively small for V and D being small; and also these are constant. So 
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in these, consistent with the classic- velocity limit v'^/c'^ ^ in question, V^V can to good 
approximation be replaced by its computed value: V^V = V[(Vy0)e*'3 + ^^^e^'^i'S/Q] ~ 
i{K + kajV^ + = —{K"^ + A;^)^, where in going to the second last expression we 
dropped the cross-term products between the mutually orthogonal Ve^^-^ and Ve*'^'^^ 
whose contribution to the final expectation value is in general zero (for an explicit proof 
see [3a,b]). Using the identity relation 7^ = 1-1-7^^, the above rewrites V^^ = — 7^/^^^. 
(iii) In the two terms ^f^V^V and c^V^V in (11), the coefficients {d'^/dT)/^ oc -iu 
and are large, with u being the scale of the particle's total energy. So, in these the V^y 
ought to be kept in functional form. But in the first of the two terms the large d^/dT 
itself effectively will be unaffected by the small V and D, and can therefore be replaced 
by its computed value —iuo'^ (used the small Q assumption), thus = —i2Duj. 

Substituting with the reduced forms of (i)-(iii) for the respective V^^, and 
in (11), simplifying using the basic relation i^^7^c^ = o;^ and the relation mc^ = hw 
given in (A.l), multiplying the resulting equation by — ^, (11) finally reduces to 

ih^ = V'y + + i2L»W2y + i/JiD^^-^V. (16) 



m " 1^1 



We next proceed to separate in wave equation (16) the K- and the K^- processes, 
which are inexplicitly contained in a factor 7 in each term as we will see explicitly below, 
and based on this we further reduce the equation at the classic- velocity hmit. To this 
end, with ^ formally given in (15), wc first compute each derivative in (16) explicitly, 
and expand the 7 factor (7 = l-|-i^-|-|^-|-...) in each: 



2 ' 8 c'' 



^ = = -^[^2 + f?,(l + -_ + ...)]y, where f?d = -r^d, = (-) 

Iv^u ^ _7!^ = = - ^ri + 2^ 

m 7M M ^ M 2M^ 4c2 



12, 



Vy =i(X + 7i^rf)y, |Vy|2 - (Vy*)(Vy) = {K + ^Kdf\^f. (17) 
On equal footing as the above, ^ expands in its exponent as 

y = (^gi[(i^+7^d)^-(^2+iTd(l+|^+...))r+Q]_ ^^g^ 

The condition v'^/c^ — > in general ensures K >> K^, >>> i?^- So, on the scales of 
Kd and i?^, the harmonic functions e'^^ and e~'^^ oscillate so rapidly that they present 
to any external observation effectively constants. Hence, e~*^^ ~ 1, e*^^ ~ 1; and 

lim y = = Z^Z^B = IZ^, = Ce^I^'^^-^'^^l. (19) 

Taking accordingly the classic- velocity limit of the relations of (17), substituting in the 
resulting relations with (19) for ]P and its derivatives (V^ip' = —K^'I^, ^ = —i ff^ 
= iKd^, V^* = -iKd^*, iVlZ^p = i^Jl^Z^P for small Q assumption as earlier) for 
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the K^-, f^d- terms while keeping the K-,f2-teTm.s as computed values which are large 
and will be unaffected for V and D being assumed to be relatively small, we have 

lim = -if2W+ — , lim = ti- + ^t, lim = iK^+V^, 

lim VT ^ -iK^ + ^^\ lim I P = X^I^^P + I V^Z^P. (20) 

We dropped the cross-term products in the last relation of (20) for similar consideration 
as earlier. Finally, subjecting wave equation (16) to the classic-velocity limit and 
substituting in the resulting equation with the expressions of (20) we have 

hQ^ + ih— = -—^ - -TT W + Vjlf - i2DhK'^^ + iDhV'^^ 
dT M 2M 

+!^^ [i^vi-^i' + i-^ri v-z^i^] w. (21) 

Equation (21) multiplied by ^ contains a component equation 

hf] = - i2DhK^ + ipiDhK^ (22) 

for a monochromatic electromagnetic wave produced by the given source but at zero 
velocity, and is not of our direct interest here. This equation holds always true for a 
given particle of a fixed rest mass and can be subtracted from ^x(21); multiplying ^ 
back to the resulting equation from left, we obtain 

ih— = -^vV + Vc^ + iDhV^^ + ipiDh^—-^^. (23) 

If V varies arbitrarily with X, thus V — V{X,T), cp^ and </?* are in general no 
longer plane waves. On the other hand, assuming V{X, T) is well behaved, we can 
divide L into a large, number of small divisions of width AX each. In each small 
division, {Xj,Xj + AX), the potential, V{Xj,T) = Vcj, continues to be approximately 
constant and is exactly so in the limit AX = 0, and here the above plane wave method 
holds valid. Elsewhere, V{Xj,T) = 0. Going through therefore the foregoing procedure 
similarly for each division, j, with j = 1, . . . , X, we obtain equations of identical forms 
as (19), (23), etc., except with K^, f^d etc. denoted by ^Kaj, Kdj, ^dK^p etc. The 
{^Kdj(^iT)ys are mutually orthogonal and form a complete set. So the total wave 
function is the sum 

T) = ^Y. T) = ^Y. A^^Ce-^""'^^-^'^ ■ e^^^^^; (24) 

or, 

^X 

with AK^.Ce-'^<ii'^+'^ = 27rXlf=i '^{Xs,T)e-'^''y^' the Fourier transform of !Z^(X„T). 



- / AT ' ^ 



^N 
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Multiplying -^Aj^^. through the corresponding equation of (23) for ^j, summing 
the equations over all j values we have 

d^^A^ = -7^7 E A.,.^. + E E A.,.^. 

dT 2M ^ * ^ ^ \fN ^ * 

j j j 

+^^^v ^ E + ^/^^^^ (^)^EAV. ^ 

(25) 

Where, = ■ ■ ■ + • F(X,_i, T) + 1 ■ V{X,,T) + ■ l^(X,+i, T) + . . . = l^(X„ T); 

A^^. = Aii^j since the amplitude of the physical displacement ^ must be real; p — 
E,- E,- ^L/i = E,- E/ for and mutually orthogonal and A*,,, 

real; and Ej Ej I VA^,.!Z^,f = V J^M^dj^^^^ '^'^dj^j mutually 
orthogonal. 

Substituting (24) in (25) wc obtain a generalized result of (23), a wave equation 
describing the kinetic motion of the particle in an arbitrarily varying, well-behaved 
potential V: 

ih— = 7T W + V^ + iDhV^^ + il3iDn}- (26) 

dT 2M p ^ ^ 

Equation (26) is seen to represent an ordinary Schrodinger equation except for the extra, 
nonlinear term iDhV'^^ + iPiDh^-^^^ due directly to the frictional force /. 

5. Diffusion currents. Continuity equation. Doebner-Goldin Equation 

Making some standard algebra to equation (26) and its complex counterpart leads to 
an equation for the total current jtot = jqm + Jobs of the probability density p = jip'p: 

^ + V(jq^ + jobs) + iPi - l)2D^-^\^\' = 0. (27) 
Where 



Jqm = ^iii^^*)"^ - "^^^n Jobs = -DVp = + 



(28) 



with jqm the usual quantum diffusion current and jobs the observable diffusion current 
as earlier except now expressed in terms of the classic- velocity limit function ^. The 
first quantity, jqm, has an imaginary diffusion constant Dqm = ^ and this we know is 
to an external observer non-observable. 

Suppose there are no "sinks" in the medium nor external reservoir in contact 
to the medium that trap or conduct the total current jtot- So the particle and the 
(continuous) medium as a whole is adiabatic — a condition having an equal footing with 
the "unitary representation of vector field (of diffeomorphisms group)" employed in [2]. 
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Then, equation (27) for the total probabihty density current of the particle, of a wave 
function ^ governed by wave equation (26), needs to conform to the continuity equation 
which, for D being real and jobs being observable, is of the Fokker-Planck type: 

^ + V(jq„+iobs) = 0. (29) 

Comparison of this with (29) suggests the third term in (27) must vanish; so /3i = 1. 

With the A value in (7) we find: W + W* = ^f^, W = +2^(V!^*)!Z>, and 
W* = +i^iy^)^* . With the /3i value in turn directly in wave equation (26), we 
finally have 

or 

th— = H'^, H'^H + iDhG, H^-—V' + V, G = VV+ (^^).(30a) 

Equation (30) is seen to be exactly the Doebner-Goldin form of nonlinear Schrodinger 
equation, the Doebner-Goldin equation, introduced in [2]. In view of their respective 
meanings, the unitary representation of vector fields in [2] and the probability density 
conservation in an adiabatic total system here are apparently two alternative but 
equivalent conditions. It is thus natural that the use of the latter here has led to 
the same result as based on the former in [2] . 

As was well appreciated in [2] (1994), the nonhnear total Hamiltonian H' as of 
(30a) is in general complex and not Hermitian as would be required by the usual 
linear Schrodinger equation; a complex nonhermitian Hamiltonian is today a topic of 
increasingly many studies. In this regard the foregoing derivation of equation (30) based 
on the lED particle model additionally points to that, underlining the complex H' and 
its imaginary part iDhG respectively are a complex total force -F^icd ~ -^mcd + ^|/mcd| 
and an imaginary frictional force /mcd, a property drawn the author's attention by D. 
Schuch at the SNMP conference, Kiev, 2007. As is suggested by the mathematical form, 
we may comprehend the imaginary /med as a physical variable orthogonal to the real 
Fined- As such, a measurement of the total force F^ed would then inform the modulus of 
it, iF^edl = V-^med + /med' ^"^^ direct addition of two scalar component forces 
and also not an ordinary vector sum 'P'^^^^ = Fmcd + fmcd- 

Concerning the solution for the Doebner-Goldin equation (30) we shall later only 
refer to an interesting and also relevant case treated by H.-D. Doebner and G.A. 
Goldin in [2]. Starting with the denotations specified in (24)' we have the more 
general expressions: p = l^pe^^Q,^ ^ ^ _2\^\2^-^Q2d^^ ^^^^ = ^[(Ve)^ + TV^ + 
2(Vg2)|eP]e-'«^ Jq^ = 2fe[(Vr)C - - 2^(Vgi)|eP]e-2«^ Following[2] we put 
Q2 — 0; the foregoing then become: 

pHer. 1-1^; (31) 

j,» = ^l(vr){-rv{-2,(voi)ijp], = A|(v{-)5 + rva. (32) 
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Secondly, suppose the system is in stationary state, of which one of two possible 
descriptions is ^ = 0. (Since for Q2 — 0, p of (31) does not contain D explicitly, 
so the "stationary state" here is not a small-D approximation but is exact as long as 
(30) holds. But we derived (30) based on a small D condition, which agrees with the 
small D requirement in [2].) Then p = of (31) is independent of time. Combining 
this with (27) follows V(jqm + jobs) ~ 0. Or, jq^ = DVp + B with B a constant. 
Substituting in this last equation with (32) for jqm, restricting V[(V^*)^ — ^*(V^)] = B 
to be independent of time which ensures that when = 0, ^ is a solution to the 
usual Schrodinger equation, and further with the specific choice S = 0, one gets 
2^(— 2ipV(5i) = —DVp. Or, VQi = —^^^ where F = mD/h. Integrating gives: 
Qi — — r'ln l^p. Substituting in (24)' with this solution for Qi and the Q2 — earlier, 
one gets: 

^^Z%, ^^Ce-'"^"^, 2; = e-'^i°l^l'. (33) 

Paper [2] also discussed that the other of the two possible descriptions of the stationary 
state to be = 0. This will also find a significant application in the examples later. 

Different forms of the nonlinear term would imply other boundary conditions than 
an adiabatic one, or other applied forces than of the form here. In recent years, in 
terms of group theoretical approach H.-D. Doebner and G.A. Goldin [2] (1994) and 
A.G. Nikitin and A.G. Nikitin and R.O. Popovych [6] gave classifications of nonlinear 
Schrodinger equations in association with diffeomorphism group representations and in 
general terms. D. Schuch gave an insightful review [7] on the nonlinear Schrodinger 
equations proposed by different authors with analysis regarding the quantum physical 
justifiability of solutions, and introduced an interesting logarithmic form of nonlinear 
Schrodinger equation. H.-D. Doebner, A. Kopp and R. Zhdanov generalized in [6] 
nonlinearity to Dirac systems. Corresponding representations of these and beyond based 
on the lED particle model in the future can be similarly of value for test of the model and 
for gaining insight into the corresponding mechanical nature of nonlinearity of quantum 
systems. 

6. Damping in dielectric media as generic application of the 
Doebner-Goldin equation 

In most applications the motion of a macroscopic object will in general be dissipated 
or more restrictively, damped, to a greater and lesser degree. The dissipation is 
typically known in the form of heat exchange with the environment and manifesting 
as an observable diffusion current. But such a description for a single quantum particle 
system, as described by (30) being in stationary state, needs be taken in an effective, 
average way only. This directly follows from the circumstances that heat reflects in 
general an energy current composed of many random collisions of a large population 
of individual (quantum) particles, during which the particles in general deviate from 
stationary state. Apart from its possible applications in the aforesaid effective way, it 
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has been desirable[2, 4] to know whether there may exist a Doebner-Goldin form of 
observable diffusion accompanying a stationary- state quantum particle literally and as 
an intrinsic process. The lED particle model underlining the foregoing derivation of 
the Doebner-Goldin equation implies in fact the presence of such processes applicable 
essentially to all quantum particles in any dielectric media; we elucidate these below. 

Consider an lED particle is moving in the total medium of an ordinary material 
medium n and the penetrating vacuum taken here literally to be dielectric, of a total 
dielectric constant k as measured against a true empty space, the space after removal of 
the dielectric vacuum. An explicit knowledge of the structure of the dielectric vacuum^ 
is not needed for the dielectric relations given in this paper. The total medium and the 
particle are as a whole evidently adiabatic. Measured against the true empty space, the 
particle's component radiation electric field propagated in the total dielectric medium is 
E, and would be E^ if "propagated" [8] in the empty space. (In this section we shall for 
conciseness drop the superscript j, either because this is not directly of concern or the 
variables actually may represent the geometric mean quantities.) When measured in the 
usual way against the vacuum with the vacuum regarded as effectively "non dielectric", 
the field E in the total medium would be E and represent the same force (note 
that the charge involved apparently causes no effect, for an original discussion see [31 
e, g]) acting on the same medium as measured in the same inertial frame and must 
therefore be equal, E = E^; this is irrespective of against which medium the force is 
measured and represented. Supposing for simplicity the material medium is nonpolar, 
with the vacuum being naturally nonpolar, so the total dielectric medium is nonpolar 
and the charge produces in it a depolarization field Ep. The corresponding dimensionless 
wave displacements accordingly are: y(= = E/A, = E^/A, and (P = E^/A. 

Applying the standard dielectric theory for ordinary materials to the generalized 
dielectric system of an ordinary material and the vacuum here we can write down 
the following relations (for a systematic treatment see [3g,e]): E(= E^ — ^) = 
Ep^-XE,E^E, + Ep, with 

k = kX^ X = K-l = {x'l + l){x, + l)-l. (34) 

Where, x is the susceptibility of the total dielectric medium and Xo that of the pure 
dielectric vacuum each measured against the empty space; is the dielectric constant 
and Xn the susceptibility of the ordinary material medium n measured in the usual way 
against a "non-dielectric" vacuum; e^^, (= Cq) = /toe0, is the permittivity of vacuum and 
the permittivity of the empty space. Multiplying by I /A, taking the classic- velocity 
limit as in (18), the dielectric relations for the fields in the above become then 

■J'^-X^, + (35) 

I There exist today various propositions for the contents and structure of the vacuum as held in different 

fields like in QED, QCD, etc., or by individual authors inchiding the "vacuuonic vacuum structure" 
proposed by the present author [3g-h,e]; there appears to exist no direct experimental information 
regarding the explicit structure of the vacuum. 
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Comparing the first relation in the above with the general form of wave function for the 
Doebner-Goldin equation, W — Z^^ of (19) or more generally (24)', we have 

Z = 1/k. (36) 

This states that, the damping factor Z corresponds rather generally to the inverse of 
the dielectric constant k of the medium in which the particle resides. In the case where 
p is independent of time, as specified by the Q2 = and small D conditions, ^ and 
T are described by the specific solutions (33), which combined with (36) gives the 
corresponding expressions for the two dielectric parameters 

«: = l/2; = e^^i^l^l', x = e^^'"l«l'-l. (37) 

In the specific case when no ordinary material presents, wc have a pure dielectric 
vacuum, thus = 1, x° = 0; (34) and (35) reduce to k, = Kq, x = Xo — '^0 " 
Wo = ^ =W^ + 3^0 = -Xo^o; and (36) and (37) reduce to Zq = 1/kq. Zq and 
are for a specified particle here evidently universal constants, given that the vacuum is 
ubiquitous, isotropic and uniform throughout the space to as far as we know all of the 
time. As a consequence, the wave function W^ appears to have never directly manifested 
itself in our present day's detections which are commonly based on the variation of 
the wave amplitude of a particle as a function of location and time; our only direct 
knowledge of the particle wave appears to be the Wq (= Wq ) of which the Zq or Kq is an 
inseparable component. Despite this, wc sec that first of all there presents a complete 
agreement between the prediction from the Doebner-Goldin equation, applied to the 
lED particle, that the electromagnetic waves "inside" (or comprising) a particle can 
in general admit damping but without decaying with time in amplitude, and the fact 
that electromagnetic waves "outside" (i.e. detached from) a particle, becoming directly 
observable, factually essentially do not decay with time in amplitude in the vacuum. 

Further, the dielectric vacuum, hence the Ej^ field of a charge in it and accordingly 
the W^ wave in the empty space, has an indirect yet profound manifestation according to 
a recent theoretical prediction [3j,f] by the author with coauthors. Namely, the Ep field 
participates to produce an attractive depolarization radiation force acted universally 
between two particles 1,2 of masses Mi and M2 and charges qi and q2, separated at a 
distance R. This force is as the result of the Lorentz force in their mutual E^, -Bo(= -B°) 
fields say in the case of a pure vacuum: Fjj/ = qj^^^-^jf^, — 1,2. The geometric 
mean of the mutual forces is Eg — a/] < E12 >< E21 > \ = '^^V^^ , where <> represents 
time average, — e, C — irxoe^ /^Ih^Pi, e is the elementary charge and the other 

constants are as specified earlier; this force Eg was elucidated in [3j,f] to resemble in all 
respects Newton's universal gravity. To this application of the Ep field, the present study 
adds that the Epi^ field of particle i producing the depolarization radiation force leads 
directly to a damping in the particle's Schrddinger wave iZ^^j, by a factor Zq, and the 
associated extra Hamiltonian term is a Doebner-Goldin nonlinear term added to that of 
the usual Schrodinger equation. In Sec. 7 we shall explicitly express the force directly 
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responsible for the damping. Also in this context, the jqm = solution mentioned earlier 
may be a case where the particle wave and thus jqm is shielded, say by a material wall. 
And on the other side of the wall j^m = 0; here only the Doebner-Goldin observable 
diffusion current jim prevails. This directly corresponds to the property of the gravity 
which can not be shielded by any materials and on the other side of the wall as here it 
will propagate alone. 

In another specific case when an ordinary dielectric medium presents and we 
represent the vacuum in the usual way as non-dielectric which thus effectively plays 
the role of an empty space in the dielectric relations, the total dielectric medium thus 
reduces to the ordinary dielectric material medium alone, thus Kq = 1, Xq = 0. The 
relations of (34)-(35) now reduce to k = k^, k° — 1 = ip'n = etc., and (36) 

reduces to 1/kI = Finally, with = (37) reduces to < = 1/Z° = e*^^" 

Xn = e*^''^'^"'^ — 1. That is, represents now a damping of the wave function IPq 
one would measure in a pure vacuum medium, into one will measure in the ordinary 
material medium of a dielectric constant k.^. This is an usual representation of a particle 
in a material medium; this is directly comparable to the familiar fact that a radiated 
(detached) electromagnetic wave in an optical material in general experiences a complex 
dielectric constant and susceptibility. 

The above two specific situations would in general simultaneously enter in our 
material world, where a material particle commonly is more or less surrounded by 
other material particles or substances and which, in the extreme case when all ordinary 
material substances are absent, is left to be the dielectric vacuum itself. Besides, a 
material particle is electromagnetic in nature (in the sense that such a particle invariably 
contains an electric charge), which is a direct observational fact and is irrespective 
of the specific lED particle model employed in this study (whereas the lED particle 
model only is essential in leading to a formal relationship between the wave function, 
the electromagnetic field of the charge and the corresponding depolarization field of 
the particle). These two universality features of the material world determine that 
a depolarization (radiation) field presents intrinsically universally with a particle. It 
therefore follows that a Doebner-Goldin damping, identified here with a depolarization 
radiation field, is an intrinsic phenomenon presenting always to a material particle. Such 
a prospect that the nonlinear process could be intrinsically universally accompanied with 
the particle process as elucidated in the above two examples was strikingly conjectured 
in the Doebner-Goldin original paper[2]. 

7. Self depolEirization radiation force: Gravity from the dielectric medium 

We now give a concrete expression for the frictional force due to a total dielectric 
medium, of dielectric constant k, against a particle moving in it. In the medium the 
particle's component radiation fields are , B^{= E^ /c), j = ] for the fields propagated 
in the direction parallel with v and j — % for fields antiparallel with v similarly as earlier. 
These fields are the results after damping in the medium, from the un-damped El, 
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measured in empty space, by a depolarization radiation field — —{El — E^) — —xE^ 
and a corresponding = —{Bl — B^) = —x^B^ due to the presence of the medium, 
where Xm — ~ 1) ^^^^ c = Cb/a/k (the permeability is assumed to be 1 here). The 
charge, due to the fields induced by itself, is acted by a magnetic force according 

to the Lorentz force law: 



where = qE^/M; sign= + for j — ] and = — for j — %■ Equation (38) expresses that, 
irrespective of the sign of the charge and of the momentary directions of the alternating 
fields generated by the charge to its right (i5^,±i?^ with a velocity c) and to its left 
(-E*,=F-B* with a velocity — c), F^^ is always a pull to the charge on either side from 
the medium. In this connection, F^.p refers to a self depolarization radiation force on 
the particle; and this represents a gravitational force on the particle from the dielectric 
medium. 

While the particle is in stationary state, its (oscillatory) charge is constantly 
traveling in an alternating -\-X- and —X- directions. The charge when traveling in 
the —X- direction is similarly acted by a pull on either side, but now the "sign" = + 
for j — X (wave generated opposite to the charge motion direction) and "sign" = — 
for j = f. So, on average the charge is acted from either side by a pulling force, an 
attraction, given by the geometric mean of the two Doppler-displaced forces: 



Fm.p = ^JfLpFLp = (sign) ^^^^'^' x, (39) 

where E — y/E^E^. This mapped to the medium is similarly a pull to the deformed 
segment in question by its surrounding in the medium, corresponding to a reduced 
displacement of the medium, this is opposite to the tensile force associated in general 
with the usual E^ field. 

We are here mainly interested in the resistance produced by Fm.p against the 
total motion of the particle, that is the time rate of F^.p as measured over a certain 
time interval AT: J^'^ {dF^.p/ dT)dT ~ ATdF^p/dT. With this, the frictional force 
following the usual definition is: 



ATrfF^ p AT ,dF^._„ dFrr^ ^dX 



Jm.p - ^ - + QX Ot' 



where E = and W = dX/dT as earlier. Similarly as Fm.p, fm.p is always opposite in 
direction to the tensile force T^. We see that indeed, in both its acting as a resistance 
against the particle total motion and in its functional form, f^.p of (40) resembles fully 
the frictional force / expressed formally in (8) earlier. 
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Abstract. We set up the Maxwell's equations and subsequently the classical wave 
equations for the electromagnetic waves which together with their generating source, an 
oscillatory charge of zero rest mass in general travelling, make up a particle travelling 
similarly as the source at velocity v in the field of an external scalar and vector 
potentials. The direct solutions in constant external field are Doppler-displaced plane 
waves propagating at the velocity of light c; at the de Broglie wavelength scale and 
expressed in terms of the dynamically equivalent and appropriate geometric mean 
wave variables, these render as functions identical to the space-time functions of a 
corresponding Dirac spinor, and in turn to de Broglie phase waves previously obtained 
from explicit superposition. For two spin-half particles of a common set of space- 
time functions constrained with antisymmetric spin functions as follows the Pauli 
principle for same charges and as separately indirectly induced based on experiment 
for opposite charges, the complete wave functions are identical to the Dirac spinor. 
The back-substitution of the so explicitly determined complete wave functions in 
the corresponding classical wave equations of the two particles, subjected further 
to reductions appropriate for the stationary-state particle motion and to rotation 
invariance when in three dimensions, give a Dirac equation set; the procedure and 
conclusion are directly extendible to arbitrarily varying potentials by use of the Furious 
theorem and to three dimensions by virtue of the characteristics of de Broglie particle 
motion. Through the derivation of the Dirac equation, the study hopes to lend 
insight into the connections between the Dirac wave functions and the electrodynamic 
components of simple particles under the government by the well established basic laws 
of electrodynamics. 
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1. Introduction 

P.A.M. Dirac established in [1] a relativistic quantum mechanical wave equation, 
Dirac equation, for a point electron based on the relativistic energy-momentum 
relation subjected to Lorentz transformation under rotation. In [1] P.A.M. Dirac also 
theoretically predicted for the electron the existence of an internal oscillation state, 
a magnetic moment, and by interpretation of the negative energy solution, an anti- 
particle state known today as the positron. The Dirac equation has proven to be an 
accurate equation of motion for (two) spin-half quantum particles at high velocities; 
most notably, Dirac predicted based on his equation the relativistic intensities of Zeeman 
components of spectral lines and the frequency differences [1 (1928b)] in exact agreement 
with experiment. Up to the present however it has remained an open question that 
what is waving with the Dirac wave functions, or Dirac spinor, a similar question as 
for the Schrodinger wave functions and the de Broglie waves ? In addition, the Dirac 
theory meets with a few its own open questions. What is the nature of a Dirac internal 
oscillation? How are the Dirac space-time functions explicitly connected with the spin 
orientations, the signs of charges, the signs of the energies, and in the extreme situation 
when an electron and position annihilate, the emitted two gamma rays and conversely? 
What is the symmetry of the total spin of an electron and positron? Also, in the case of 
an isolated single electron or positron in zero external field where the spin orientation 
is of no consequence, it would be desirable to have a way to directly write down the 
corresponding Dirac equation without involving the Pauli matrices. These as well as 
various other not fully addressed questions relating to fundamental physics seem to 
consistently point to the inadequacy of the point particle picture of today and the need 
for a representation of the internal processes of the particles. 

Recently, using overall experimental observations as input data the author proposed 
an internally electrodynamic (lED) particle model [2a] (with coauthor P.-I. Johansson) 
or sometimes termed a basic particle formation (BPF) scheme, which states that a 
simple (basic) particle like an electron and positron, etc.. briefly, is constituted of an 
oscillatory point-like (elementary) charge with a specified sign and a zero rest mass, 
and the resulting electromagnetic waves in the vacuum. As a broad test of the lED 
particle model and also as an endeavour of understanding the various puzzles relating to 
fundamental physics, in terms of solutions for the electrodynamic processes of the model 
particle with its charge's sign and total energy as two sole input data, the author has 
further achieved with coauthor(s) derivations/predictions of a range of basic properties 
and relations of the simple particles [2a-j] including the relativistic mass, de Broglie wave, 
de Broglie relations, Schrodinger equation, Einstein energy-mass relation, Newton's law 
of gravity and Doebner-Goldin equation, among others. As to the Schrodinger wave 
function specifically relevant here, the solution[2a,c] showed that it is the (envelope of 
the) standing wave, superposed from the Doppler-differentiated electromagnetic waves 
generated by the particle's travelling source charge, that is waving. 

As previously shown e.g. in [2c], the direct solutions for the classical wave equations. 
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derivable from the Maxwell's equations, for the electromagnetic waves comprising a free 
particle consist of Doppler-displaced plane waves; these superpose to two opposite- 
travelling beat waves that resemble directly the de Broglie phase waves and in turn 
the Dirac space-time functions which in common are functions of the particle's total 
energy and linear momentum and thus are " relativistic" . It therefore is foreseeable 
that the classical wave equations for the electromagnetic waves would more naturally 
lead to a wave equation of the particle corresponding directly to the Dirac equation in 
comparison to the Schrodinger equation. We elucidate in this paper a formal procedure 
which transforms the classical wave equations for the electromagnetic waves of two spin- 
half particles, of identical space-time functions and tending to approach one another, to 
the Dirac equation. Through the procedure we show that the Dirac internal oscillation 
corresponds to the oscillation of the electromagnetic waves at a geometric mean of 
frequencies which in general are Doppler-displaced owing to source motion, and we 
elucidate the explicit relationships between the internal electromagnetic waves, charges, 
spins, the centre-of-mass and total wave motions and the associated energies of the 
particles under the government of a few established elementary laws of electrodynamics. 

2. Wave equations for the electromagnetic waves of particle. Solutions 

We consider an lED particle, here an electron or positron, is as its source charge q {= e 
or — e) travelling at a velocity v in +2;-direction for the present along a one-dimensional 
box of side L in the vacuum. The charge q of the particle has an oscillation associated 
with a total energy Sq, which is minimum at f = 0, denoted by 6^; 8.q may be endowed 
e.g. in a pair production in the vacuum. In virtue that it describes the ground state, 
8,q cannot be dissipated or detached from the charge except in a pair annihilation. 

The charge q of the particle generates owing to its oscillation electromagnetic 
waves of radiation electric fields E-^'s and magnetic fields B-^'s described in zero applied 
potential field by the Maxwell's equations as: 

V • E^' = Pl/eo, V • B^' = 0, V X B^' = ^of, + {l/c')dt^\ V x E^' = -dtB^. (41) 

Where p^^ is the density and j;^ the current of the particle's charge, assuming no other 
charges and currents present; eg is the permittivity and /iq the permeability of the 
vacuum, and c is the velocity of light; dt = d/dt. Expressing the jth fields generally 
by a dimensionless displacement ip\ = Dip\ thus = E^ /c = Dip^/c, with D a 
conversion constant, considering regions sufficiently away from the source only so that 
= = 0, and with some otherwise standard algebra of the Maxwell's equations (41), 
we obtain the corresponding classical wave equations for the electromagnetic waves (/t^ 's 

c^vV^' = ^tV^ (42) 

with = d'^/dt^. In the above, j — ] labels the component wave generated in the 
direction parallel with -\-v, and j — X the wave parallel with —v] within walls there 
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prevail also their reflected components generated by the reflected charge at an earlier 
time and being at the present time as if generated by a virtual charge travelling in the 
— z-direction, labelled by j = virf and j = vir|. j is to distinguish a Doppler effect 
owing to the source motion to be expressed in (47) below. In Appendix A we outline in 
relevance to the particle model a few further standard relations of classical and quantum 
electrodynamics for the electromagnetic waves, and a derivation of the particle's mass 
given by the author previously[2a,e] (with P.-I. Johansson). 

To the particle we now apply an electromagnetic force F = Fg + F^, with 
Fg = — ?V0a the Coulomb force in z- direction and F^ = — gv x V x the Lorentz 
force due to an external scalar potential (pa and vector potential A^, expressed in SI 
units as for all other quantities in this paper. F^ may be simplified using the BAC-CAB 
rule as F^ = — g[V(v ■ A^) — Aa{v ■ V)]. In the applications below A^ is constant in L 
or in each small division in question (see end of Sec. 3), and v is constant in L for the 
particle being in stationary state and also is parallel with V and z, so V(i' ■ A^) = 
and F^ = qAa{v ■ V) = qAavV. Thus, F = -qV4>a + qAavV. The formula of F is in 
the usual usage estabhshed for a point particle; so when extending to the extensive lED 
particle here, F apparently directly acts on the particle's point charge. 

We need to map the F to a force directly interacting with the internal fields , B^, 
or (f^ of the particle. We observe that, in virtue of its form, (42) represents just a classical 
wave equation for a mechanical wave of a transverse displacement aip^ propagated in an 
apparent elastic medium, a being a conversion factor of length dimension and apparently 
being cancelled in (42). On grounds of this direct correspondence, but taken as a 
heuristic means only in this paper (so that we here need not involve the details of this 
elastic medium), F therefore interacts with the internal fields through a force F^^^ 
directly acting on this apparent medium. We can think of the medium to be composed 
of coupled dipole charges which do not move along the 2;-axis but the F^^^ propagates 
across the dipoles at the wave speed c. If viewing in a frame where F^^^ is at rest, then 
effectively the dipole charges are travelling at the speed c; so as a first step of mapping, 
the Lorentz force on the medium ought to scale as F-''^ = {±c/v)Fm = itgAacV with 
+, — for the j = t, I waves; thus F-'' = Fg + F^'m- Under the actions of the respective 
forces, the acceleration F^' /m^ of the particle's charge of a dynamical mass (due 
to the charge's total motion and equivalently the ip^ motion, see further Appendix 
A), and that of the medium of a dynamical mass 971^, -^med/^^j must equal, i.e. 
F^'/m^ = FUm^- Thus F^, = ^F^' = ;^(-V0. ± A.cV). 

By its pure mechanical virtue the force F;^^^ acting on the continuous medium is 
nonlocal and will be transmitted uniformly across the medium here along the z-axis 
of effective lengths L^^,L^^ for the j = |, J waves (</?•' winds loops about L). Using 
the geometric mean = a/L^LJ,, thus V0a = ±(0a/L^)i and V = ±1/L^. With 
these, putting — p^LP^ where p^ is the (geometric mean) linear mass density of the 
medium, writing for conciseness V0a and also the final F^^^j in scalar forms and keeping 
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the generally arbitrarily oriented Aa in vector form only, F^^^ becomes 

^med = -P<^V^'' V = q<Pa-qAaC, V* = - gA„c. (43) 

^med ^® implemented in (42) by directly establishing the corresponding wave 
equation for the apparent elastic medium acted by -F^ed- without F^^^^, the elastic 
medium would be deformed owing to the disturbance of the oscillation of the source 
charge alone, by a total displacement u = (i^,j^{,t. and be thus subject to a tensile 
force Fr = P;C^. The applied F^^^j and add up to a total force acting on the particle 
through acting directly on the medium 

F'n -Fn- FL^ = p, [c^ + V^/m^] . (44) 

Where, the minus sign of F^^^ is because this force tends to contract the chain. 
Assuming ip^ is relatively small which in general is the case in practical applications, 
F^' is thus uniform across the L. A segment AL of the medium along the box, of 
mass A3Jt^ = /SSJl^/ P — Pi/^L ~ Pi^z, will upon deformation be tilted from its 
equilibrium position z-axis an angle -d^ and + Ai?-' at z and z + Az. The transverse 
{y-) component force acting on A9Jl<^ is AF-^'^^ = F-''^[sin(i?^ + Ai?-') — sini?''] = 

F^'j^'^{a(p^)/S.z — api + ^ V^f/^^Az. Newton's second law for the mass ASD?^ 

writes p^Azdflaip^) = AF^'^^. The two last equations give the equations of motion, 
on dividing ap^Az, for per unit length per unit linear mass density of the medium at 
z or equivalently the classical wave equations for the electromagnetic waves (p^^s in the 
fields of the applied potentials 4>a, A^ 



[c' + q(<l>a - A„c)/mt] VV^ = 9,V^ [c' - q(<l>a + A„c)/m*] = d^^K (45) 

This for 0^ = = reduces to (42) given directly from the Maxwell's equations earlier. 

Assuming for the present (pa, A^ are constant and also A^ is small such that the 
particle motion effectively deviates not from the linear path, so the solution of (45) 
consists of plane waves (p^ — C/^ and (/?* = C/* (Figure la, sohd and dotted curves) 
generated in +z- and —z- directions and initially also travelling in these directions at 
speed cu^ — c, with 

yt ^ (j^i[klz-ujh+ao]^ jt ^ _(J^i[-klz+uj*t-ao] j-^gj 

(Figure 1 a-b, single-dot-dashed and triple-dot-dashed curves), G = e*^^, and C 
(= ACi/y/Z) a normahsation constant. Where, 

= K/{1 - v/c) = -f^K, k^ = K/{1 + v/c) = -f^K and = -f^f2, = ff? (47) 

are the source-motion resultant Doppler-displaced wavevectors and angular frequencies; 
~ i-v/c ^ ^* ~ i+v/c '■> — -^c are values of k^, at f = 0. (47) further gives 

kl^ = k^ -K = -f^Ka, kl = K-k^ = 'y^Kd where = {v/c) K. (48) 
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Supposing V « c (yet v^/ may be large so that dynamically the 7 factor in (49) below 
can be different from 1) and accordingly the de Broglie wavelength (A^ = 27r /{'jKd)) 
later will be much greater than the electromagnetic wavelength (yl = so at the scale 
of Xd the rapid variation of C is to an external observer no different from the constant 
1, that is lim„«c C = 1. Thus <^^|e=i = /\ ¥'*le=i = /* the /^ /* as given in (46) 
represent external-effective space-time functions. We see that P and /* are two new 
plane waves travelling each to the right at equal phase velocities, ^ — — W — — 

which is c/v times the velocity of light c. It can be checked (Appendix B) that the exact 
solutions (/7^ 's and in turn the external effective /^'s given above placed in the respective 
wave equations (45) above and (52) below yield exactly the expected relativistic energy- 
momentum relation. 

The Doppler- displaced variables A;^'s,u;-^'s in the /-^'Sjf^'s are not single valued and 
thus are not good dynamical variables of the particle. The respective geometric means 



kd^ \lk\kl=-fKd, u^V^'^iO, with 7 = = VVl - ^Vc', (49) 

are evidently good dynamical variables of particle and also are appropriate in view 
of the stochastic virtue of the electromagnetic waves. These are also the natural 
independent variables of the superposed wave functions (for a detailed elucidation see 
[2c,d]): -0 = = e^gi[fed2-(^t+ao] and V''''' = (^"'•t+(^">'t = _evirgi[fcd^+a;t+ao] ^hich are 

two beat waves of a wavelength = 27r/A;d, traveUing each at the phase speed W — /v 
to the right and to the left as due to the actual and reflected (virtual) charges travelling 
to the right and left respectively; Qd = 2Ci(^^^^+-c^''^'->^\ = 2Cie^[(-^+^'='*)^+^^*] 
and Gd=Qd"=2Cie'^^^ = 2Ci for v << c. An inspection will show that clearly ^, 
resemble directly the dc Broglie phase waves of the particle in the constant 0^, fields 
here, A^^ and kd arc the dc Broglie wavelength and wavevector, and accordingly hkd = Pv 
the linear momentum. If (pa, A^ are arbitrarily varying in L and well behaved, we can 
divide L in a large number of small divisions in each of which the plane waves remain 
true and their sum gives according to the Fourier theorem the total wave, and in three 
dimensions the de Broglie particle motion is a straightforward extension of a locally 
one-dimensional motion; the wave equations to be given will formally be otherwise the 
same (for a formal treatment in the case of a Schrodinger system sec [2c,k]). We shall 
thus for simplicity proceed the remainder of treatment for constant (pa, A^ and, until 
the discussion regarding spin rotation in Sec. 5, for a particle motion in one-dimension. 



3. Wave equation for total motion of particle 

For single particle or for many particles without regarding the spins, the functions ip 
and 0'^"', or equivalently the fr and fc of (55) below, are seen to be identical to the usual 
solutions to the Dirac equation, c.f. Appendix C. So their wave equations, originally 
the (45), evidently must have a direct correspondence with the Dirac equation. The 
remainder of the task mainly will be to identify a physically justifiable procedure to 
transform (45) to a form of the Dirac equation under corresponding considerations. 
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Figure 1. (a) shows an lED model electron constituted of an oscillatory charge — e 
of zero rest mass (g), travelling at velocity v, and the resulting Doppler-differentiated 
electromagnetic waves ip'' and 93* (solid and dotted curves) of a angular frequencies 
uj\lu^ generated in +z- and — z- directions, plotted for the real parts in a time 
interval (3/4)t in a one-dimensional box of side L — A^; r — 2Tr/uj, uj = \^uj'^uj^, 
Xd = 2tt{{v/c)uj) being the de Broglie wavelength. and /* (single-dot-dashcd and 
triple-dot-dashed curves) are the corresponding external-effective waves, shown in both 
(a) and (b). fr and fi (solid and dashed curves) in (b) are the dynamically equivalent 
mean-variable wave functions; these resemble directly the (opposite-travelling) de 
Broglie phase waves and are equivalent to the space-time functions of Dirac spinor. 



First, similarly as Dirac (or as alternatively but compatibly argued in [2c]) we want 
the eventual wave functions of the particle, and thus immediately the or the original 
if^ to be linear in hw and thus dtf-' here, so that at any initial time determines its 
value at any future time; and we want similarly for the linear momentum here. We 
shall thus transform the second order differential equations (45) to first order ones and 
in the end take the limit for c >> t; as follows. For the c^V'^if^ terms of (45), starting 
with the full wave functions (f^ = Gf^ we first lower the spatial derivative one order 
as VV^ = V[il^Kep] = i-f^K[{Ve)p + CV/^]. Restricting to v « c, we thus can 
replace VC by its computed value iKC and in turn put C=l for each term; this gives 

VV'|e=i = [-I'K^QP + i^'KQVf%=i = -^^K^P + t^^KVf, 3 = f, t- (50) 

Next, assuming (^a, relatively small as typically is true in applications, so the resulting 
force constant (i.e. force per unit displacement) on the particle does not vary across 
L; we can thus replace the V^(p^ in the (pa, terms of (45) by its computed value as 
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VV^|e=i = -i^'^K^P, thus 
^^^^VV^|e=i = q{-<f>a + A„c)2^, ^^^^^VV*|e=i = + A„c)^ (51) 

For the final expression we used Kc = Q a,s earher, and — ^^M and Mc^ = hf2 
given after (A.l) and (49). Finally, the d^ip^^s of (45) lower one order as d^ip^^lQ^-^ — 
—i'~f''f2dtf\ dfip^lQ^^ = i'~f^f2dtf^. Substituting these and equations (50)-(51) in wave 
equations (45), multiplying the first resulting equation by and the second by 

with cK = Q and UK = Mc as before and after (A.l), we eventually obtain the wave 
equations for the electromagnetic waves of the particle expressed by f \ /*: 

[Mc^ + q(j>,-c(ihV + qA,)]p = iHBtf, [-Mc2 + #„ + c(iW + gA„)]f = iHBtf. (52) 

For the particle dynamics in question we want to further transform (52) to be 
expressed by the particle wave variables, i.e. the and cu defined in (49), and the 
corresponding wave functions, the fr, fe to obtain below. We shall below obtain such 
functions through a dynamic equivalence transformation directly from the /^ /*; these 
ought to be and will show to be functions identical to the ip, ip™ obtained in a physically 
more transparent way earlier; the present approach below will advantageously preserve 
a direct tractable connection with the original /^/*, thus also v'^v^*, whose wave 
equations (52) or (45) give the relativistic energy-momentum relation exactly based on 
the Doppler equations (47), see (B.2) of Appendix B. What (accordingly) is in question 
in the transformation mainly is to maintain an equivalence to the quadratic equation 
(B.2); this corresponds to the equations |fr|^ = -§^-q^, etc., with j,j' = j, |, 
H, fx' = r,l, z" = t, z {u = 0,3); and ffM-gf = 0. The equivalence condition requires 
in particular the transformed quadratic to he = k^, that is, it has a plus sign in 

front and its cross-term product with Mc^ (i.e. the discussed after 57) is absent. This 
can be achieved if we introduce a wavevector being the imaginary of (thus orthogonal 
to) kd- 

kd = {l hl')k',. kd=i7/tf)k'„ thus fctfct = (i/,2)^2 (53) 

(compare k^ with the operator p^^p — f V later). (53) alternatively can be expressed by 

(a) : = k^k^ or (b) : {-k^M = k,k,. (54) 

We now first transform the Dopplcr-differentiated /■'(; A;^, uj^)^s (as short hand notations 
of /■'(z, t; /c^, ci;-')'s) to a pair of mean (wave)-variable functions /^(; fc^i, co')'s (denoting 
ffj_{z,t;kd,ujys) by, say, satisfying (a) of (54) and ordinarily cu — ■\/ufhj^ of (49): 

fi; kl u') ^ /.(; ka. uj) = Cr^^^^'-^'+'^°\ f (; ki a;*) ^ /,(; uj) = C,e^[^'^^+'^*+"°l ; (55) 

see these functions plotted in Figure lb. The transformed fr-, fe indeed are desirably 
identical functions to the original /^ /* if disregarding the high-order differences in 
the coefficients 7\ 7* and 7 in the wave variables and the reversed travel direction 
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of fi from /*. The fr, fe, being identical functions to the ipjip^^^ earher, indeed are 
therefore the pertinent space-time functions of the particle; these are each functions 
of the source motion and the total (electromagnetic) wave oscillation and accordingly 
directly resemble the de Broglie phase waves; and these are equivalent to Dirac's space- 
time functions. 

To entail that in the matrix representation later (Sec. 5) a cross-term product 
discussed after (57) is similarly absent, for transformation of the first derivatives we 
are compelled to satisfy the alternative condition (b) of (54). The use of (54b) and 
ordinarily the u — y/cJkJ^ of (49) first directly leads to the intermediate transformations 
for the /•''s given in the left column below: 

/^-^^-/^T^-/- ^f--iKf^^-iU-fr)-Vfr (56) 

where, // = Qe^[-^rf^-'^*-"o], /; = Cre'^-~^'i^+'^^+<\ a'^ = -ckq; the transformations (al) 
and (a2) in (56) finally naturally lead to the same fe,fr as in (55), which indeed also 
represent the original /* in all aspects (having the same phase velocities and wave 
forms) except the opposite rotating phases on the complex plane and the opposite signs 
of Oq and ao that altogether arc dynamically inconsequential. The relations in the left 
column and the use of (54b) again then lead to the results in the right column which is 
actually in question in respect to dynamical equivalence here. 

Substituting in wave equations (52) the transformation relations (55) and (56) gives 

{Mc^+q(l>a)fr+c(ihV-qA,)fi = iHBtfr, (-Mc''+qcl)a)fe+c(ihV+qA,)fr = ihdtfe. (57) 

As a check, placing in (57) the fr, ft of (55), multiplying the first and the negative 
of the second resulting equations, dividing frfe, putting for simplicity 0^ = = 
for the problem mainly of concern here, we correctly obtain the same result as (B.2): 
M^c^ - h'^ky + = where = -kj following (53) and (49); = 0. 

4. Two-p£Lrticle system: spins, charges, and time-arrows 

Consider two spin-half particles 1,2 having identical sets of {fr, feYs tend to occupy 
the same location z or more precisely the same region in (0,L); suppose these are 
noninteracting (a finite particle-particle interaction can in principle be included in 
and will not affect the general conclusions below). In virtue of the statistical 
nature of the electromagnetic displacements, the probability of finding a portion of 
particles 1 and 2 at locations Zi and Z2 is proportional to the product ffx{z\,t)f^j,i{z2,t). 
Since these have identical space-time function sets, the corresponding total space- 
time function is evidently symmetric, thus fs, here in the only form fs{zi,Z2) = 
■:^[fr{zi,t)fi{z2,t) + fi{z2,t)fr{zi,t)] to be compatible with the antisymmetric total 
spin function later. 
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In the case of two identical electrons (Figure 2, panel a), their spins then need 
according to Pauli principle be opposite (left graph in the figure) to avoid both particles 
occupying the same quantum state; the total spin function for this is antisymmetric 
(Xa = "^['^(l)/5(2) — ci'(2)/3'(l)]). Apparently it is in general also relevant that we 
introduce charge functions Q_(l), Q_(2) to reflect the sings of charges 1,2; these for the 
two like charges are trivial identities and lead to a trivial symmetric total charge function, 
Qs- The above functions together define an antisymmetric two-electron function (Figure 
2a, right graph), ip^ — fs^sXa, yielding as expected a total probability independent of 
how we sample the two stationary-state identical, indistinguishable (as result of being 
identically extensively distributed in L) particles. 

In the case of an electron and positron (Figure 2, panel b), the total two-particle 
function is symmetric (right graph), thus i/jg, as follows from the observational fact that 
two such particles can approach each other arbitrarily close and, in an extreme case 
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Q^(2)=— Q.(l) is explicitly contained 
in \|/in)'s in the graph, not in the text. 




Figure 2. Two spin-half lED particles 1 and 2 described by identical, thus mutually 

symmetric sets of Doppler-displaced space-time functions {(p\(p^}'s or effectively 
{fri fiY^ tend to occupy the same location z; panel (a) shows two electrons and (b) an 
electron and positron. Particle 1 has a spin S in +z- direction, denoted spin up, a(l), 
and is parallel with the generating direction of the fr wave; particle 2 has spin S in — 2;- 
direction, denoted spin down /3(2), and is parallel with ff, their opposite counterparts 
(termed virtual spins) /3 (1) and a (2) are parallel with fe and fr- Right graphs show 
the complete wave functions and {■i/'+(2), V". (2)} of particles 1 and 2. 

In panel (a), ipvi'i-) and ipvi^) {1^ = +,-) are antisymmetric due to the antisymmetric 
spin functions and same charges. In (b). these are symmetric (as shown in the graph) 
due to antisymmetric spins and also opposite charges (Q+(2) = — Q_(l)), the latter of 
which leads to the radiation electric fields E-' 's are opposite in direction (in the text the 
Q_^(2), Q_ (1) are not expHcitly regarded, rendering the V'y(l) and ■i/'i'(2) for the electron 
and positron to be antisymmetric and the same as for two electrons). 
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annihilate into "one point" in the vacuum. The total charge function for their two 
opposite signed charges evidently needs be antisymmetric, thus Qq. (If imagine the 
wave displacement ^p^ in the medium is executed by a chain of dipole charges then it 
is immediately clear that the corresponding radiation electric field E\ with \E^ oc ip\ 
produced by the positive charge is reversed from that by the negative charge, see left 
graph in Figure 2b.) Placing the two known functions in ips = Qa/sXsym gives therefore 
Sym=antisymmetric and thus Xsym — Xa for the electron-positron system. 

Our particles contain each the electromagnetic waves (p\(p^, or effectively fr, fe, 
generated in the specified +z- and —z- directions which assume definite relationships 
with the spin orientations that turn out in a measurement: Given, say, particle 1 is 
measured to be spin up, a{l), with fr being parallel with it, then its other internal 
process, fe, is parallel with a virtual (indicated by a prime) spin-down state, I3'{1). 
Similarly for particle 2, being then actually spin down, thus fe parallel with f3{2), its fr 
is parallel with a virtual spin-up state, a' (2). From the foregoing antisymmetric spin 
requirement follow the relations for the spin functions: 

a (2) = -«(!), = -/3(2); with /3'(1) = -«(!), a (2) = -(3(2) (58) 

following the opposite signs as is meant by the "virtual" spins. Through (58), the virtual 
spin vector of particle 1, S'(l) (virtual spin down) and the actual spin vector of particle 
2, S(2) (actual spin down), pointed each in the —z- direction, are now each represented 
as scalar quantities with minus signs. 

If disregarding the signs of charges explicitly, the electron-electron and the electron- 
positron are two equivalent systems of identical, spin-half particles, each described by 
the space-time functions /^(; kd, uj)^s with fi = r,£, which joined together with the spin 
functions give the complete wave functions: ^d) ^) = fn{\ kd, uj)ai^{n) with u = +,-, 

n = 1,2, = a, j3. Now as a further step to conform our wave equation later to matrix 
form, we hereafter require that the ipu{n) functions are elements of a matrix of one 
column, i/). The matrix wave equation itself will entail the two desired features discussed 
after equation (57), that is, (i) the product of —kd and kd in the quadratic equation is 
positive: k^ (entailed by the situation that these in the matrix form are offdiagonal 
elements, see (62) or (C.l), and (ii) the total cross-term product = (entailed by the 
characteristics that in matrix equation the '0i^(n)'s are explicitly mutually orthogonal). 
The first of these two features which we have up to now enforced by use of kd for kd, 
should no longer be used in the matrix form to avoid a dual accounting. The space-time 
functions accordingly write fr{; kd,io) = C^e*''^'*^"'^*! , fe{; kd,oj) = Qe*'*^'^^^'^*! with kd the 
ordinary scalar quantity and related with A;^, fc^ through (49). Accordingly, 

^^+(1) = a{l)fr{;kd,u) = a{l)Cr^^^^'-^'\ ^.(2) = /3(2)/,(; A;,,a;) = /3(2)ae^[^''^-*], (59a) 
V'.(l) = 0{l)fe{;kd,u:) = 0{l)Ce^^^^'+^'\ t/jp) = a {2)fe{;kd,u;) = a'(2)Qe^['=''^+-*l.(596) 

The complete wave functions (59a)-(b) (Figure 2, right graphs) describe two identical 
particles of opposite oriented actual spins and accordingly opposite virtual spins, and 
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these are identical to the solutions (see equation C.3) for Dirac equation. We can readily 
check that placing the foregoing relations in the antisymmetric total function for two 
spin-half (like charge) particles, tpa = "^["^+(1 )''/'- (2) — ip_{l)ipj^2)] correctly leads to 
that the probability of finding two identical particles at any location 2; in L is not 
altered by interchanging the locations of the particles (the indistinguishability) . We 
can also check that the same ipa is given by the product of the separate total functions: 
i^ai^i, ^2) = /s(-2i, -22)Xo(l) 2); notice that once we specified say particle 1 is spin up and 
2 spin down, then fj.{zi,t)fe{z2,t)P'{l)a'{2) and fi{zi,t)fr{z2,t)a{l)P{2) are zero since 
these do not describe the present reality. 

Lastly, the spin-up state of particle 1, a(l), is associated with an effective 
electromagnetic wave fr travelling to the right, thus dtfr/ fr = —ii-o, while the spin-up 
state of particle 2 with fi travelling to the left, thus dtfe/fe = iuj; the latter has as if a 
reversed time arrow relative to the former. We may introduce the time arrow functions 
defined for particles 1 and 2 as T(l) = 1,T(2) = —1, such that the action of these 
on the time derivatives project the wave propagations to be both in the -|-z-direction: 

5. Dirac equation 

For two identical, spin-half particles of identical sets of space-time functions fr, fe, 
described by wave equations (57) tending to occupy the same location z, we shall now 
express the corresponding wave equations in terms of the complete wave functions of 
Sec. 4. For particle 1, we thus multiply the first equation of (57) by a{l) and the second 
by /9'(1), act T(l) in front of the time derivatives, denote its charge by gi, and get 

(Mc^ + gi0j/,a(l) + c{ifN - giA„)/,(-/5'(l)) = ih7{l)dtfrOi{l), 
{-Mc^ + gi0„)/,/?'(l) + c{ihV + giA„)/,(-a(l)) = ih7{l)dt{fe(3' {!)). (60) 

For particle 2, instead we multiply the first equation of (57) by — /3(2) and the second 
by —a' (2), act both equations by T(2), denote its charge by 52, and get 

{Mc^ + g20a)/r(-/3(2)) + c{ihV - g2AJ/,(+a'(2)) = th7{2)dtfr{-f3{2)).^ 
i-Mc" + q2(t>a)fi{-a{2)) + c{ihV + g2A„)/,(+/3(2)) = iKy{2)dtfi{-a {2)). (61) 

In the second terms in equations (60)-(61) we made the replacements a{l) —0{1) 
and I3'{1) -a(l), I3{2) -a' (2), a {2) -f3{2) based on (58), to conform to the 
transformed space-time functions earlier. 

Substituting in (60)-(61) with (59a)-(b) for the tlj^{nys and the T(n)'s expressed 
earlier, and, to form a direct contrast between the actual spin directions of the two 
particles, rc-arranging the resulting four equations in the order of spin-up states of 
particles 1 and 2 first and then spin-down states of particles 1 and 2, we finally obtain a 
set of four coupled linear first order partial differential equations governing the motions 
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of the two particles in terms of ipv{n): 

{Mc^ + — c(iW — qiAa)ip_{l) = ihdtipj^l) (particle 1, spin up) 

(Mc^ - g20a)V'+(2) + c(i W + 92 Aa)V'. (2) = ihdtijp) (particle 2, spin up) 

(— Mc^ + Q'l^a)''/'. (1) — c{ihV + qiKa)ipp) = ihdt4'_{l) (particle 1, spin down) 

(-Mc^ - q2(l)a)ilJ. (2) + c{ihV - g2A„)V^^(2) = indtip, (2) (particle 2, spin down) (62) 

Prom the discussion of Sec. 4 that the ■0,^(n)'s and accordingly also their first derivatives 
are mutually orthogonal, it follows that the linear equations (62) are equivalent to a 
matrix equation. Supposing specifically the two particles are a positron and an electron 
and therefore qi — q, q2 — —Q, the matrix form of (62) is thus 

Hopij) = ihdti}, with Hop = hMc^ + qcpa + ca(p„.op - qAa) and p^.op = -iKV (63) 

being the relativistic total Hamiltonian and linear momentum operators. Where, 




(64) 

The off-diagonal elements of the matrix a^, crzii{= 1) and az22{— —1) here correspond 
to the a{l) — 1 and a (2) — —1 earher. We see that, ^ is equivalent to a Dirac spinor, 
the 2;-component of Pauli matrices, and as a whole, equation (63) is identical to the 
Dirac equation for an electron-positron system equivalent to here. Por the present case 
rotation transformation is trivial, so cr = cr^x. 

Suppose more generally the two particles' spin angular momenta, S (= |cr)'s, are 
along an axis n executing in general a precession about the 2;-axis at a fixed angle 
(arccos(^)). For each particle being in stationary state, its S (similarly its magnetic 
moment fis{— —eS/m)) as a vector quantity when in small rotations about the z-axis 
must maintain invariant with respect to its projection on the 2;-axis, n • S = if^, 
and is Hermitian. In addition to the antisymmetric condition given by the cr^ of (64) 
above, an infinitesimal rotation transformation as such needs be unitary. A specific set 
of transformation matrices having these properties are known to be the Pauli matrices, 
ax, (Jy and of the standard expressions and as expressed in (64), a = axX+CTyy+CzZ. 
And, the unitary matrix Iz about the 2;-axis naturally extends to a unitary matrix about 
the new n- axis in three dimensions, given hj I = Ix + ly + Iz. Substituting in (63) 
with a and / for and / gives a Dirac equation of the same form, now for spins in 
arbitrary directions. 
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Appendixes lA-IC: 

Appendix A. Total energy and inertia of particle wave 

As a general result of classical electrodynamics based on solution to the Maxwell's 
equations combined with Lorentz force law, an electromagnetic wave j transmits at 
the speed of light c a wave energy and a linear momentum p' — jc. Here, the 
amplitudes of e-', accordingly of p', E\ and etc., are continuous values. Following 
M. Planck's discovery of quantum theory in 1901, it has been additionally understood 
that these quantities are by nature quantized in amplitudes; an electromagnetic wave 
of frequency a;/27r has an energy e — nfuo, consisting in general of n momentum- 
space quanta, or photons, each of an energy Huj; and the classical continuous amplitude 
solutions to these are only approximations when n is large. In the present problem, 
in conformity with experiments, especially the pair processes, the electromagnetic wave 
comprising our basic particle has a "single energy quantum", n = 1; so e = hu. It has 
been further proven especially through quantum electrodynamics that the Maxwell's 
equations, and the subsequent classical wave equation (42) or (45), continue to hold, 
and the quantisation of the fields and wave energy etc. is the result of subjecting 
the canonical displacement and momentum, the u{= a(p) and ii here, to the quantum 
commutation relation [u, ii] — ih. 

The total wave of our particle of a single "quantum energy level" fuj in a one- 
dimensional box has, following the solution to the Maxwell's equations earher [see after 
(42)], two components, (p^ and </?*, with their frequencies being Doppler-displaced to cu^ 
and cj* as a result of the source motion as given in (47), which are related to u through 
(49). For the total wave comprising the particle, e represents therefore a dynamical 
variable of the particle, here the total energy of the particle. 

The electromagnetic waves, E\B^^s or v^-^'s, rapidly oscillating at frequencies 
(c<j-'/27r)'s, of a geometric mean frequency lj/2tt and wavelength A = c/ {u/2'n-) will, when 
ignoring the detailed oscillation as will effectively manifest at some distance, appear as if 
being two rigid objects, wavetrains, travelling at the speed of light c. In view that their 
speed of travel, c, is finite as contrasted to infinite, the wavetrains have inevitably each 
finite inertial masses, m-^'s, thus an inertial mass m = \fm)m} for the total wavetrain 
and hence its resulting particle. This mechanical depiction of the total wave, as a 
rigid "wavetrain", permits us at once to express according to Newtonian mechanics the 
linear momentum of the wavetrain to be p = mc. Combining this with the classical 
electrodynamic result e = pc above gives the kinetic energy of the wavetrain e = mc?, 
being equivalent to the Einstein's mass-energy relation. This energy and the Planck 
energy earlier ought to equal, thus 

m = huj/c^; or at t; = 0: M = h(2/c^ (A.l) 

with M the rest mass of the particle; combining (A.l) with (49) gives m = 7M. 
Combining (A.l) with p = mc further gives mc = {hLj/c^)c = hk and accordingly 
Mc = hK, with a; = /cc, i? = Kc and k = jK as earlier. 
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Appendix B. Relativistic energy— momentum relation for the 
electromagnetic waves of particle 

Consider first the simpler case of Aa — 0. Placing in wave equations (52) with f\ /* of 
(46), dividing the resulting first and second equations by and — /* and sorting give 
Mc^ + hk\c — hu;^ — q(f)a, Mc^ — hk\c — fuj^ + g0a. Multiplying gives 

- h^k^ky + Q = h^uj'uj' - q^(t>l {B.l) 

Where klk^ = kj and u^u^ = u'^ following (49); Q = Mc'^ch{k\ - k\), with - = 
2A;d(^)7 ^iid -^c^ — ^Kc, so Q — 2h^kd'c^. With these, putting hk^ — ±p„, hu; — ±e 
where py,e are here variables having positive and negative solutions and thus the right 
hand side of (B.l) reduces as ■\/[{huJ — q(t)a){hjj + q(j)a)Y = a/[— (^ — q<t>a){^ — Q'^a)]^ = 
(e — g0a)^, then (B.l) reduces exactly to M^c^ + p^c^ — ~ Q'PaY- This, or this in the 
more familiar form for 0^ = 0, 

gives just the experimentally widely corroborated relativistic energy- momentum 
relation. For the more general case of finite, denoting /c^ = k\ — fc^' = k\ + 
the particular feature that (the cfi^cctivc portion of) A^ is always perpendicular to kd^P 
leads to kf = k'Jkl' = kj- q^Klin\ or, (ip'J^ = {±hk'^f = Tihk^ - qAa){-hka - 
qAa) = [t{Pv - ^Aa)]^. (B.2) thus generalises to M^c^ + c^Pv'^ = - Q(l>a)'^- 

Appendix C. Solution of Dirac equation from the standpoint of particle 
internal process 

We shall here mainly discuss the choice of the solution forms of the Dirac equation 
from the standpoint of internal processes of the lED particle model for simplicity for 
spins along z-axis, in an otherwise basically standard procedure. The two equations 
of (60) or (61) for particle n = 1 or 2 are coupled in ipj^n) and ip_{n) and can not 
be solved separately as in Appendix B. We need to solve each two, or more generally 
the four equations of the Dirac equation (60) together. Let the trial functions be: 
■0^(72,) — Csne^^'"^~^*\ s — n — 1,2. Placing these in (63) and rearranging give 



/ £-Mc2-#„ -[p„-gA„]c \ 

£ - Mc^ - q(f)a [p^ - qAa]c 

-[p„-gA„]c e + Mc'^-q(Pa 

\^ [p„ - qAa]c e + Mc^ - q(f)a J 



( n^) \ 
V'-(i) 

V V'.(2) J 



(c.i) 



(C.l) corresponds to the four linear, homogeneous algebraic equations (C.2) for the 
ips,/s as four unknowns below; for these to have nontrivial solutions, the determinant 
for the matrix of the coefficients of (C.l) needs be zero. This is det= [{e — q4>af ~ 
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M^c^]^ — c^p'^ = 0, with p'^ = Pv — ^Aq. This has two degenerate sets of square roots 
solutions: e — q<j)a = ±1/ M^c^ +pu'^c^, which being identical to (B.2). In view that each 
particle has internal processes, we thus naturally assign symmetrically two of the four 
solutions to particle 1, as e — q(j)a = ±a/ f\/Pc^ + Pv'^c^, and the other two for particle 2 
as e — q4>a = T\/ M'^d^ + pj'^c^- These two distinct sets of square-roots solutions to the 
algebraic equation above represent two (distinct, identical) particles, like an electron 
and a positron, which do not transit from one to the other, a point agreeing with 
reality and having been stressed by P.A.M. Dirac from the very beginning in [1]. These 
algebraic solutions are in contrast to the usual problem of eigen values arising generally 
from boundary conditions and being each possible states of same particle between which 
transitions generally can occur. 

With and e as known parameters, we further solve the four algebraic equations 

(£_(Mc^ + ^0j)^^(l)^p;cV'.(l) (a), {e-{M^ + q4^a)W)^-M_{'^) (b) 
{e+{Mc'-q<l>a))^_{l)^p',c^Jil) (c), {e+{Mc'-q<l>,))^_{2) ^ -p'^c^p) {d){C.2) 

corresponding to (C.l), for the wave functions. Taking the imaginary of equations (a) 
and (b) first, multiplying the resulting equation with equation (c) and the second with 
(d) respectively on opposite sides, substituting with tp*{l)ip_{l) = Cl^, i/'*(l)i/'_^(r 
V'*(2)V'.(2) = C\ and V;(2)V'+(2) = Cj, we get 
{e-{M^ + q(t>a))Cl = (£ + Mc2-#„)C_2i, - (Mc^ + #„))C| = (s + Mc^ - g^jC^^ 
These have two independent solutions, and in mathematical terms two of the four wave 
functions can thus be arbitrarily chosen. In view of the lED particle model by which each 
particle has internal, wave processes consisting of two components in the one- dimensional 
box, it is natural here that we choose the values for and symmetrically, in the 
sense also = — C^, with these in the two equations above the values for C_i and C_2 
then follow to be uniquely given as 



-( — 



e+{Mc^-q(f>. 



'a)) 



1/2 



c.l 

C_2 



( / £-(Mc2+#a)" V^^ ^ 



where |C_|iC_i| = C^, |C^C_2| = C^. With the above in the trial functions, we get the 
complete solution for Dirac equation 

1/2 



( H^) \ 

V ^-(2) ) 



e+{Mc^-q4>a) 
e-(Mc^+q(I>a) 

e+{Al(:--q6a) 
e-(Mc^+q<pa) 



e-(Mc2+q4>a) \ 
e+{Mc^-q4'a) J 



(J ^i{kdz-ujt) 



\ 



1/2 



1/2 



\ 



:-{Mc^+q(i>a) 
:+(Mc2-q(ha) 



1/2 



(C.3) 



/ 



Agreeing with the wave functions directly based on lED particle model, ^_^(1), (1) are 
two opposite travelling component waves of particle 1, and ■iIjJ2),iIj_{2) of particle 2; in 
the meantime, the spin-up component waves of particles 1 and 2, ipj^l) and ip_p) travel 
in opposite directions and similarly the spin-down component waves. 
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